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Abstract: 

We compute the mixed correlation function in a way which involves only the 
orthogonal polynomials with degrees close to n, (in some sense like the Christoffel 
Darboux theorem for non- mixed correlation functions). We also derive new 
representations for the differential systems satisfied by the biorthogonal polynomials, 
and we find new formulae for the spectral curve. In particular we prove the 
conjecture of M. Bertola, claiming that the spectral curve is the same curve which 

appears in the loop equations. 
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1 Introduction 

Consider a pair of random hermitean matrices Mi and M2, of size n, with the proba- 
bility measure: 

This random 2-matrix model has many applications to physics (in particular in quan- 
tum gravity, i.e. statistical physics on a random surface and conformal field theory 
[^1^) and mathematics (bi-orthogonal polynomials 05). Another important appli- 
cation of the 2-hermitean matrix model comes from the fact that it is the analytical 
continuation of the complex-matrix model, which describes the Dyson gaz aX (3 = 2 
and is an important model of Laplacian growth jHH]. And the complex matrix model 
plays a crucial role in the AdS/CFT correspondance, in the so-called BMN limit |31j . 
In that model, it is important to know how to compute mixed expectation values. 

For all applications, one would like to be able to compute various expectation 
values. Some expectation values can be written in terms of eigenvalues of Mi and M2, 
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for instance < Tr Tr M2 >, which we call non-mixed because each trace contains 
only one type of matrix Mi or M2 but not both. On the contrary, mixed expectation 
values are those where Mi and M2 may appear together in the same trace, for instance 
< Tr M^M2 >. Mixed expectation values cannot be written in terms of eigenvalues of 
Ml and M2, and are thus more difficult to compute than non- mixed ones. 

Beyond the technical challenge of computing them, mixed expectation values should 
play an important role in applications to boundary conformal field theory and to the 
BMN limit of string theory/CFT correspondance. While many formulae for non mixed 
expectation values have been known for a long time, in particular in terms of bi- 
orthogonal polynomials [ini EH 1211 121 HI HI U] ; it is only recently that formulae have 
been found for mixed traces. In particular the following expectation value: 



was first computed in jH]. The idea was to diagonalize the hermitean matrices Mi = 
VXV^ and M2 = VUYU'^V^ where U and V are unitary matrices and X and Y are 
diagonal matrices containing the eigenvalues of Mi and M2. Then: 



and using the Morozov's formula [331 HH] for unitary integrals of the form < UijUli >, 
one reexpresses Wn in terms of eigenvalues of Mi and M2 only. Then the integration 
over eigenvalues is done with the method of biorthogonal polynomials jUl- The result 
found in jS] is thus a. nxn determinant involving recursion coefficients of bi-orthogonal 
polynomials: 



where the matrices Q and P implement the recursion relation (multiplication by x and 
y) of the biorthogonal polynomials, and n„_i is the projection on the polynomials of 
degree < n — 1. Those notations are explained in more details in section |231 

The mere existence of such a formula was a progress, but a n x n determinant is 
not convenient for practical computations and for taking large n limits. In the non- 
mixed case, the Christoffel-Darboux theorem allows one to rewrite expectation values 
in terms of smaller determinants, whose size does not grow with n 01 Ej . 

The purpose of the present article is to find a similar property for Wn{x,y), i.e. 
write it in terms of determinants whose size is independent of n. 

As a byproduct of such a rewriting, we are able to find new formulae for the spectral 
curve of the isomonodromic differential system satisfied by the corresponding biorthog- 
onal polynomials. 




(1-2) 




(1-3) 




(1-4) 
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1.1 Plan 

• In section 2, we recall the definitions of bi-orthogonal polynomials, as well as their 
Fourier-Laplace transforms, Cauchy transforms, and the third type of solution 
introduced by |T3l 1^ . We also recall their recursion relations, obtained by 
multiplication by x or derivation. We define the infinite matrices (finite band) Q 
and P, which encode the recursion relations {Q and P are also Lax matrices). 

• In section 3, we study the inverse oi x — Q and y — P, in particular, we show 
that they have distinct right and left inverses. The difference between the right 
and left inverse, is related to the so-called folding matrix [TT], and is found to be 
a sum of bi-orthogonal polynomials and their various transforms. 

• In section 4, we introduce the kernels |l6^, which are the building blocks of corre- 
lation functions. We show that they satisfy some Christoffel Darboux relations. 

• In section 5, we compute the biorthogonal polynomials, their Cauchy transforms, 
and the kernels as determinants involving matrices Q and P. 

• In section 6, we introduce the notion of windows of consecutive bi-orthogonal 
polynomials jHl EH E] , because the recursion relations (matrix Q and P) allow 
to rewrite any bi-orthogonal polynomial pm as a linear combination of pj with 
n — d2 < j < n, where d2 + l is the degree of the potential entering the weight of 
orthogonality. 

• In section 7, we prove our main result. We compute the mixed correlation function 
Wn{x, y) in terms of polynomials in the window, and in terms of the kernels. 

• In section 8, we study some consequences of that formula. In particular, we 
find new representations of the differential system satisfied by a window, and we 
compute the spectral curve. We find several new formulae for the spectral curve, 
and we prove the conjecture of M. Bertola, which claims that the spectral curve 
is the same curve which was found from loop equations. 

• In section 9, we discuss the consequences of that formula in terms of tau functions. 

• In section 10, we give some explicit examples of our formulae, namely the gaussian 
and gaussian elliptical cases. 

• Most of the technical proofs are put in the appendix. 



3 



1.2 Main results 



This article is very technical, and the purpose is to give many formulae for effective 
computations with bi-orthogonal poljTiomials. We propose some new formulae in al- 
most every paragraph. Let us mention here the most important ones: 

• In section we give integral representations of the right {R) and left (L) 
inverses of x — Q. One of the important relations is Eq. 

- Rix) = ^oo(x) ^U^) (1-5) 

• In section we give determinantal expressions of bi-orthogonal polynomials and 
their transforms, as well as kernels. For instance the kernel Kn{x,y) is: 

Kn+iix, y) cc det ({x - Q){y - P')) (1-6) 

nxn \ / 

• In section [7[ we state one of the main results of this paper, i.e. some formulae 
for the 2-point mixed correlation function Wn{x^y): 

1 1 



Wn{x,y) = l + (Tr 



X- Ml y- M2 

det 1 



nxn 
,,2 



X - Qy - P 



7n^n+i{x,y) Jn~i{x,y) X 

det (1 - nr^^i - M)f/.(i - 

(1-7) 

and theorem 17.21 gives Wn in terms of kernels only. 

We also find interesting recursion relations for Wn+\ — Wn- 

In section ISl we compute new representations of the differential system Vn{x) = 
'^'^{x)'^n{x)~^ , and we compute the spectral curve: 

£nix,y) = idetiy - P„(x)) = -7^ det(l - t/^f/*) = det(y - H^'^\x,x)) (1-8) 

The most interesting is that we prove Bertola's conjecture j3]: 

£nix,y) = iV;ix)+y)iV^iy) + x)-n+(TT- ^ 



X- Ml y- M2 

(1-9) 

This conjecture has important consequences in terms of tau functions. Indeed the 
Miwa-Jimbo-Ueno approach of isomonodromic tau functions [2Zll2HlEni, general- 
ized in , allows to express the tau function in terms of residues of the spectral 
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curve, and this formula is particularly convenient for that purpose. It shows 
that the tau function is the matrix integral, and it shows that some additional 
parameters could be added to the model. 

2 Definitions and notations about bi-orthogonal 
polynomials 

This section recalls well known facts about bi-orthogonal polynomials, and stands 
here just for setting notations and describing known properties. Notations are similar 
(although with small differences) to those of [TUl lllj. 

2.1 Measure and integration paths 

Consider the weight: 

= e-(^i(")+^2(^)+"J'), (2-1) 

where Vi is a complex polynomial of degree di + 1 and V2 a complex polynomial of 
degree d2 + I- 

di+l d2+l 

V^{x) = J2ikx' , r2(y) = ^4/. (2-2) 

fc=0 k=0 

We write the leading coefficients of V( and V^: 

t={d^ + l)ta,+i , t=(d2 + 1)4+1- (2-3) 

We choose a basis of di contours 7*^*^ with i = 1, . . . , di, going from 00 to 00 in sectors 
where the integral J e~^^^^^dx is convergent, and we choose a basis of ^2 contours 7^*^ 
with i = 1, . . . , (i2, going from 00 to 00 in sectors where the integral / e~^^^'^^dy is 
convergent (see [IHIE])- 

Then we choose a dual basis of di contours 7*^*-* with i = 1, . . . , di, going from 00 to 
00 in sectors where the integral J e~^^^^^^dx is convergent, and we choose a dual basis 
of d2 contours 7 with i = 1, . . . ,d2, going from cxd to 00 in sectors where the integral 
J Q+V2(y)(iy jg convergent, such that: 

^»n7(^-) = 5,^. , 7»nf''^ = 5,,. (2-4) 

Then, we choose did2 numbers Kij such that at least one of them is non- vanishing, 
and we define a path F: 

r:=EE'^..V^^X7(^'^ (2-5) 
1=1 j=i 
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and we define the following measure on F: 



(2-6) 



Remark 2.1 Generalized path T and matrix models 

We have introduced the generalized integration contours F, because it is the most general 
contour on which the measure dfj, can be integrated. It corresponds to a generalization of 
the hermitean 2-matrix model. Indeed, hermitean matrices have their eigenvalues on the real 
axis, and the hermitean 2-matrix model Eq. corresponds to the case F = R x R . 

A generalized path F can also correspond to a matrix model, with matrices which are not 
hermitean. It corresponds to normal matrices (i.e. which can be diagonalized by a unitary 
transformation, but with complex eigenvalues), with pairs of eigenvalues constrained to be 
on F. This allows to define an ensemble of matrices which is noted Hn x Hn{T), see |18j for 
more details. 

In this normal matrix model, it makes sense to compute matrix expectation values, in 
particular the mixed correlation function Wn{x, y) = 1 + (Tt -^ttm^ y-M2 



More definitions 
We define (see [IS1I2S1): 



ho= c//i(x,?/) = y V/ti W dx / rfy e-(^i(^)+^^ 



'2{y)+xy) 



And: 



-g^Hy) ■- 



2i7r J^(i) 



2m 



e^^'e^i^^) dx 



i = 1,. . .,d2 



i = 1, . . . ,di 



g^'^ := 0, 

^(°) := 0, 



which are the independent solutions of the differential equations: 

V^{d/dx)g'^'\x) = -xg^'\x) 

v;{d/dy)g^Hy) = -yg^'\y)- 

Then we define the following "concomitents" I25j: 



—— / dx dx 



2171 j^o) 



X — X 



W— dy dy 

2l7r J;y{i) feO) 



y-y' 



)y 



QV2{y')-V2iy) ^{y'-y)x_ 



(2-7) 

(2-8) 
(2-9) 

(2-10) 
(2-11) 

(2-12) 
(2-13) 



The c*^*-^'''s (resp. c'^^-'^'s) are independent of y (resp. x). Due to the dual choice of 
contours Eq. (j2-4|) . they are normalized: 



(2-14) 
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Indeed, integrating by parts, we can replace both V({x) by y and V({x') by y. If 
i ^ j, the contours 7^*^ and ^^^^ do not intersect, and the integration by parts gives 
no boundary term and the result vanishes. If i = j, the two contours intersect and we 
have a boundary term. A way to compute it, is to write the pole — x') as the 
sum of a principal part and 2m5{x — x'). The principal part is integrated by parts and 
gives zero as in the i ^ j case, whilst the (5-term corresponds to the boundary term in 
the integration by parts, and it gives 1. 

This computation was first done by the authors of fSl 12^1; in the comparison of 
the two Riemann-Hilbert problems derived for biorthogonal polynomials pTlj ESI I2nj 

2.2 Bi- orthogonal polynomials 

The monic bi-orthogonal polynomials [SBIISI], (if they exist), are uniquely determined 
by: 

Pn{x)=x^ + ... , g„(y)=y" + ..., (2-15) 

and 

Pn{x)qm{y) dl2{x, y) = K ^nm- (2-16) 



r 

For given potentials Vi and V2, bi-orthogonal polynomials exist for almost every choice 
of r (in fact they don't exist only for an enumerable set of F's, see |inj). 

2.2.1 Wave functions 

We define: 

^„(a;):=^p.(x)e~^^(^) , ^^[y) ■.= ^q^[y) e'^^iv) . (2-17) 



'"n 



2.2.2 Cauchy transforms 

We introduce the Cauchy transforms P ITU]: 

^„(y) := ^e^^(^) / ^p„(x')rf/x(x',y'), (2-18) 
V rin Jr y ~ y 

0„(x):=^e^^(^) / -^q^{y')dfi{x',y'). (2-19) 

2.2.3 Fourier-Laplace transforms 

We also introduce the following functions [16] IT7j: 

^^\y)-= [ Mx)e~"'dx fori = l,...,di, (2-20) 
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<P'i\x):= (t>n{y)e-^''dy for^ = l,...,(i2. (2-21) 
They are the Fourier-Laplace transforms of ipn and (j)n- For i = we will also write: 

(2/) ■■= My) , 0f (a^) := (2-22) 
2.2.4 Third-type functions 

The authors of ^31 llEj have introduced the following functions: 

^»(x) := ^ ^ rfy / d^{x\y') ^^^lM^^e^^(^)e^^p.(x'). (2-24) 
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We will also write: 

(x) := V^„(x) , 0f (y) := 0„(i/). (2-25) 

Formaly, those functions are the "inverse Fourier transforms" of the V^'s, as de- 
scribed in [101123 ■ 

2.3 Semi-infinite vectors and matrices 

We introduce semi-infinite vector notations: 

V'oo(a;) = (^o(a;),V^i(a;),z/'2(a;),...)* , 0oo(z/) = (0o(z/), 0i(2/), 02(z/), • • •)*> 

(2-26) 

and more generaly: 

i^J\x) = (4^)(x),#(x),#(x),...)* , ^JHy) = i<p^o\y), <Pi\y),^2\y), ■■■)'■ 

(2-27) 

And 

^oo{y) = {i)Q{y),i)i{y),'tp2{y), ■ ■ - f , 0oo(x) = (0o(a;),0i(a;),02(a;), . . .)*, 

(2-28) 

and more generaly: 

i^J\y) = {i^f{y).i^f{y)Af{y)....r , ^^J\x) = (^^^^(x), 0«(x), 0«(x), . . .)*• 

(2-29) 

We also introduce the basis vectors: 

n oo 



(0,... ,0,1,0,...)*, (2-30) 



i.e. the vector whose only non vanishing component is in the n + 1 position. It is such 
that: 

el^ij^{x) = iJn{x). (2-31) 
Similarly, we consider the projection matrix: 

n+l oo „ 

n„ = diag(C?^,0) = X^e.e*, (2-32) 

j=0 

with n + l ones on the diagonal. It is the projector on the span of eo, ei, . . . , e„. We 
also define: 

II" = 1 — n„_i = projector on the span of n, n -|- 1, . . . , oo, (2-33) 

and: 

m 

= Bje* = n^IIm = Ilm — n„_i = projector on the span of n, . . . , m. (2-34) 



j=n 



We also introduce the following cxd x (^2 + 1' 



(resp. oo X (di + 1)) matrices: 





:=(W°^ 






(2-35) 




:=(0j°) 




. 0j'^)), 


(2-36) 


^ oo 


- ( r 


r (1) 


■ Voc j ' 


(2-37) 


<l 


— 

•- 10OO 


0OO 


• • 0OO ) ■ 


(2-38) 



2.4 Recursion relations for the bi- orthogonal polynomials 

It is well known^ that we have the following recursion relations [T^ lUj: 

Xj^n = Qnmi^m , y4>n = -Pnm0m, (2-39) 

m 



V^n = X] Pmni^m , 'P'n = Qmn<Pm- (2-40) 

m m 

Since Vi and V2 are polynomials, Q and P must be finite-band matrices, i.e.: 

Qnn, 7^0 iffn-rf2<m<^+l, (2-41) 

Pnm 7^0 iffra-di<m<n+l. (2-42) 
In vector notations we have: 

xi)oo = Qi^oo , y4>oo = P4>oo, (2-43) 

V^;, = PVoo , 0'oo = QVoo. (2-44) 



and acting on a polynomial, gives a polynomial which can be decomposed on the basis of 
biorthogonal polynomials. 
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2.5 Relations between Q and P 

The matrices Q and P have the following properties [13 El UH] '■ 

[g,P*] = Id. (2-45) 



Qn,n+l = Pn,n+1 = \/ := 7n+l- (2-46) 

-P* + VY(Q) is a strictly lower triangular matrix 

{P' + V;m+ = , (P* + K(Q))„,„_i = - (2-47) 

7n 

and Q* -|- V2(P) is a strictly lower triangular matrix 



n 



(Q* + y2'(P))+ = , (Q* + l^2'(^))n,n-l = -• (2-48) 

7n 

2.6 Recursion relations for the Cauchy transforms 

The Cauchy transforms also satisfy recursion relations pUj : 

yifooiy) = P'^oo{y) + eo , x0oc(x) = Q'<foo{x) + eo, (2-49) 

0o(y) V^o(a;) 

^:„._g,-^ + i-SM^e.. (2-50) 

2.7 Recursion relations for the Fourier-Laplace transforms 

For z 7^ we have PUDl: 

yVoo = P ^oo , a;0oo = <5 0OO , (2-52) 



^1^) = -Q^«.^^^ , 0'« = -P0J^^ (2-53) 

Notice that they satisfy the same recursion relation as the Cauchy transforms without 
the non-homogeneous term. 
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2.8 Recursion relations for the third type functions 

For z 7^ we have UHl HSl 123 : 



- + y^,/.e--e- - Jeo, (2-54) 

^l^) = P*V^oo^^^-^?«eo, (2-55) 

y<PJ ^ = P<PJ^+ Q/Qy _ J eo,", (2-56) 

0l^) = g*0j^)-^«eo. (2-57) 

Notice that they satisfy the same recursion relation as the wave functions, with an 
additional non-homogeneous term. 



3 Inverses 

The formula for mixed correlation functions found in [8|, is written in terms of the 
inverse operators of x — Q and y — P, thus we study them in detail in this section. 
X — Q and y — P also have distinct right and left inverses, which were shown to play a 
crucial role in the notion of folding onto a window in [2] • 



3.1 Inverse 

By definition, the infinite matrix l/{x — Q) has elements: 

( — ^) = J-^ [ Qmiy') ^-^Pn{x')dfi{x',y') 

\X — ^ / nm yil'nilm JV X — X 

= 1pn{x)(i)m{x) + Rnm{x), (3-1) 

where 

Rnmix) := / qM) P-^""^ - P-^""'^ df^ix', y') (3-2) 

yilnllm Jv X — X 

which is a polynomial in x of degree n — 1. 
In vector notations we have: 

— ^ = i^oo{x)(t)J {x) + R{x). (3-3) 
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Similarly 

= 0oo(i/)^"oo*(i/) + R{y) (3-4) 

y- P 

where 

Rmn{y) := / (3-5) 
3.2 Right inverse 

The semi- infinite matrix R{x) (resp. R{y)) is polynomial in x (resp. y) and is strictly 
lower triangular: 

Rnm = if m > n (resp. Rnm = if m > n). (3-6) 
It is a right inverse [TT] of {x — Q) (resp. {y — P)): 

{x - Q)R{x) = Id (resp. {y - P)R{y) = Id ). (3-7) 



But it is not a left inverse, we have 



i?(x)(x - g) = Id - (resp. i?(y)(y - P) = Id - M^). (3.8) 



Notice that ^ 

Rn,n—1 ~ ~ Rn,n—1- (3-9) 

In 



Notice that we have: 



X — x' 



J2Rnm{x)i^m{x') (3-10) 



(resp. '^"^^^ '^^^y'^ =-J2^nUy)My'))- (3-11) 

y-y ^ 

In particular at x = x': 

^'^{x) = -R{x)^^{x) = PVoo(x), (3-12) 

(resp. (t>L{y) = -R{y)<Poo{y) = gVoo(y)). (3-13) 

Notice also that: 

R{x,) Rix,) = -^^^^1^^ (3-14) 

Xi - X2 

and 

[R{xi),R{x2)]=0. (3-15) 

It can be found, by solving directely the system (x — Q)R{x) = 1 for a lower 
triangular matrix R, that: 

det ( nr+\ (x - Q)nr+\ ) for m < n - 2 

Rnmix) = { ^ for m = n~l (3-16) 
for m > n 
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3.3 Left inverse 




Figure 1: Definition of tlie contours 7^ (y) and 7_ (y). 



Consider y G 7'-*-', tlien deform the contour 7*"*^ so tliat it crosses 7*-*''at y, and define 
7_ (y) the part of 7 wliicli stands on tlie riglit of 7*^*^ and which ends at y, and 7_,_ (y) 
the part of 7 * which stands on the left of 7*^*^ and which starts at y (see fig^). Then 
define: 



LnAx) ■■=77^J] [ dy<Pm{y) e-^^ f ^ dy" 
2t7r ^ \M\y) 



+ dy"]My")e^y" (3-17) 



It has the following properties: 

Theorem 3.1 Lnm{x) is a polynomial in x. 

L is upper triangular, with Lnm{x) = if n + d2 > m. 
L is a left inverse of x — Q: 

L{x){x-Q) = l (3-18) 

We have: 

Lnmix) = (^) + (3-19) 



proof: 

The proof is in appendix A □ 
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Alternative definition 1 

It can be found, by solving directely the system L{x){x — Q) = 1 for an upper triangular 
matrix L, that: 

det ( (X - Q)K~f^ ) for m>n + d, 
(x) = { ^ for m = n + d2 (3-20) 

for m < n + (^2 

This proves that there is only one left inverse of a; — Q which is upper triangular. 

Alternative definition 2 

Define: 

where yi{x) are the solutions of V^iy^^x = 0. Notice that /(x) is an upper triangular 
matrix such that: 

lnm{x) =0 a n + d2 > m. (3-22) 
The left inverse ^T] of x — Q is given by: 

-l{x) + l{x){V^{P') + Q)l{x) 

~l{x){Vi{P') + Q)l{x){Vi{P') + Q)l{x) + . . . . 

(3 - 23) 

Since l{x) and V"2(P*) + Q are strictly upper triangular matrices, each entry of that 
infinite sum, is actualy a finite sum. L{x) is a strictly upper triangular matrix, such 
that: 

Lnm{x) = if n + d2 > m. (3-24) 

This definition gives clearly a left inverse of x — Q, and since the upper triangular left 
inverse of x — Q is unique, it must coincide with the first definition. 

Remark: I computes the Euchdean division of qm by V^'(y) + x: 

My) = -(^2(y) + ^)Y1 ^nm{x)My) + (deg < ^2 - 1). (3-25) 

n 

Similarly, we define: 

~ ^-^ R\xi{y)) 
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and the left inverse of y — P is: 

1 



L{y) := -l{y) 



I + {Vim + P)l{y) 
= -l[y) + l{y){yi{Q') + P)l{y) 

-iivWHQ') + P)!{y){V;{Q') + P)l{y) + 

(3 - 27) 

3.4 Relationship between right and left inverses 

We have: 



x-Q , 

= L„„(a;)-5]^«(a;)0»(a;) 

i=l 



i.e. in vector notations 



1 {i)t t 



x-Q 

This imphes: 



L{x) - R{x) = J2 i^oo^\x)cfJ (x) = *oo 

i=0 



which is a matrix of rank ^2 + 1- 

Since R is lower triangular, we have, for n <m: 

Lnm{x) = J2^ni^)4>';!;!{x) if n < m, 

i=0 

and since L is upper triangular, we have ior n > m — + 1' 

Rnm{x) = - 5^ Vi')(a;)0«(a;) if n > m - ^2 + 1- 

1=0 

Notice that if n < m < n + ^2 — 1, we have: 

d2 

i)^^\x)^^^{x) = ifn<m<n + d2-l, 

i=0 

and thus: 

1 \ _| 'il^n{x)(i>m{x) \in<m 
^-Q)nm~\ -E£i4'H^)0^H^) ifn>m-d2 + l 
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4 Kernels 



The kernels are the building blocks of non-mixed expectation values |2] , and we 
need them because they will appear in mixed expectation values as well. Thus we recall 
their definitions and some properties below. 

4.1 Definition of the kernels 
4.1.1 K kernels 

We introduce the kernel: 

n— 1 

Kn{x,y) := = 0oo*(l/) n„_i ^oo(a;) (4-1) 

and more generaly: 

n-l 
m=0 

= <pj'^\y) n„_i i;J^\x) + Ki;''\x, y) (4-2) 

where 

:= (4-3) 

irj^'°)(x,y):=-— e-^^(^) / -^e^^^^"^e^'ydx' (4-4) 
2171 J^(i) X — x' 

(a;, y):=-— e'^^^^) f e^'^^'^ e'^' dy' (4-5) 

2tTC J^(^) y-y' 

K^^''\x,y) := TTT^ /" dx"j rfy" /" rf/i(x', i/')e''^(^'V^(^") e^^"+^"^ 

1 W)-W)^ 1 Vi{x')-Vi{x" 



{x — x') {x — x") y' — y" {y — y') {y — y") x' — x" 

(4-6) 



4.1.2 J kernels 

We also introduce the kernels: 

n-l 



J„(x, y) := Jo{x, y) - ^ (pm{x)il'm{y) ■= Jo{x, y) - 0oo(a;)*n„-iV'oo(y) (4-7) 



?Ti=0 



and more generaly 



n-l 



4^-^-)(x,y) := t'\x,y) - J2 (4-8) 



m=0 

16 



where 



Mx,y) Jt'\x,y) e^^^^^e^^^^) / ^ dfi{x' ,y') (4-9) 

Jy X — X y — y 

j(0'*)(x,t/) :-e^i(^) / -^e-^'^^'^e-^'^dx' (4-10) 

4'^'\x,y) ■.= c-''y (4-12) 
Sometimes one writes abusively: Jn{x,y) = Ylm=n^rn{x)ipm{y) = ^oo^n^Voo- 



4.1.3 i7 kernels 

We also introduce the kernels: 



n-l 



Hi:'\x,,X2) ■.= J2^^r^i^^)^^ri\x2)+H^''\x,,X2) (4-13) 
m=0 

where 



xa-xi 2z7r J^O) J;^(i) y-y' uj{x2,y) 

yi(.0-vi(-.) f f , ,Vi{y)-Vi{y') 1 1 a;(x',i/') 



pl/Ua:ij-\/Ua:2j /• /■ 

hS'''\xi,x2) := — 2^^^ — y^^^.^c^i/ y^dx^ 



^w. 1/ - 1/' a;' - x' - X2 Lv{x2, y) 

(4-15) 

i/^'°^(a:i,X2) :=0 (4-16) 

i/(°'°)(^i,^2) ^^e^i(^i)-^^(^^) (4-17) 

X2 - Xi 

Notice that at Xi — X2 we have: 

lim {X2 - x,)Hi''^\x,, X2) = 5ij (4-18) 

X1—*X2 

4.1 A H kernels 

Similarly: 

n-l 

Hi''\y, y') ■■= E + H^o'\y. y') (4-i9) 



m=0 
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4.2 Christoffel— Darboux matrices 

We define the Christoffel-Darboux matrices 0: 

An := [n„_i, Q] = 'jnen-iei - (1 - n„_i)Qn„_i, (4-20) 

Bn := [n„_i, P] = 7ne„-ie^ - (1 - n„_i)Pn„_i. (4-21) 

Notice that since Q (resp. P) is finite band, An (resp. i?„) is non-vanishing only in 
a sub-block of size ^2 + 1 x ^2 + 1 (resp. di + 1 x di + 1) , i.e. 

{An)ij ^0 iff n - ci2 < ^ < n and n - 1 < J < n + ^2 - 1 (4-22) 

{Bn)ij ^0 iSn-di<i<nandn-l<j<n + di-l (4-23) 
We say that An and Bn are "small matrices", i.e. their size is not growing with n. 

4.3 Christoffel-Darboux theorems 

Using the recursion relations Eq. ()2-39|) to Eq. ()2-56p . we have the following Christoffel- 
Darboux theorems [HI EI : 

Theorem 4.1 

(X2 - x,)H^:'\x,,X2) = (fJ^\x,)An^J^\x2) (4-24) 

(d^, + d,,)H^:^\xi,X2) = -(fJ^\x,)Bi^J^\x2) (4-25) 

{x - d/dy)K^:^\x, y) = <pJ'^\y)An'iljJ'\x) (4-26) 

{y - d/dx)K^:^\x,y) = <pJ^\y)B'^^J^\x) (4-27) 

^ (i)t - (j) 

{x + d/dy)Jn{x,y) = iy)AniJoo (x) (4-28) 

^ (i)t -^ (■!) 

{y + d/dx)Jn{x,y) = <j)J' {y)Bl^^J"{x) (4-29) 

In particular, at Xi = 0:2 = x in Eq. ()4-24j) . and using Eq. ()4-18|l . we recover the 
duality of ^ : 



(i)t, 



ix)An^ljJ'\x)=6,,, (4-30) 
and we find: 

Hi'^\x,x) = <fJ^\x)An^'^\x). (4-31) 
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5 Determinantal formula for the orthogonal poly- 
nomials and kernels 



5.1 Determinantal formulae 

The following formulae are very useful to express kernels or polynomials, as well as the 
Cauchy transforms, as n x n determinants. The proofs are given in appendix B. 

1. • determinant of n„_i(a; — Q)Iln-i'- 

det(n„_i(x - Q)Un-i) = Pnix) , detiUn-iiy - P)n„-i) = Qniv) (5-1) 

2. • determinant of n„_i^j7^n„_i: 

det(n„_i T^nn-i) = -j== <Pn-i{x) (5-2) 

x-Q ^/h^i 

1 Q-V2(y) . 
det(n„_i -n„_i) = ^n-i(y) (5-3) 

y-P 

3. • kernel Kn'- 

det(n„_i(x - Q){y - P')Un-i) = /i„ e^^(^)+^^(^) Kn+i{x,y) (5-4) 

4. • kernel J„: 

1 1 Q-Viix)-V2(y) 

det(n„_i — = 7 Jn-i{x, y) (5-5) 

x-Qy -P* K-i 

5.2 Inverses 

It is usefull to compute also the inverses (i.e. all minors) of the previous matrices. The 
proofs are given in appendix ^2 The following formula give inverses of n x n matrices. 

1. • inverse of n„_i(a; — Q)Iin-i'- 

{iin-i{x - Q)n„_i)-^ = n„_i - n„i?n„_i (5-6) 

2. • inverse of n„_i^n„_i: 

\ -1 " * 

n„_i^-n„_i ) = n„„i(i + ^!^^i^n")(x - g)n„_i (5-7) 



x-Q 



'-'n-l 
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3. • inverse of n„_i(a; — Q){y — P*)n„_i: 

(n„_i(a;-g)(y-p*)n„_i)-i = n„_ii?*n„ - -^^^0^*^ n„i?n,_i (5-8) 



n— 1 



4. • inverse of n„_i^^^n„_i: 



n„_i(y - p*)(n„_2 + ^(1 - n„_2)^^oo0oo*(i - n„_2))(x - g)n 

(5-9) 

5. • equivalent formula for Yin-i^;zQ:ippt^n-i'- 

""-^x-gy-p*""-^ 
_ 1 1 J 

a; - g 1/ - ^ 

= n„_i — —Yin-2 — 777n„_i + J„_in„_i'?/'oo0oo*n„_i 
2; - g y - -P 

(5 - 10) 

6 Windows 

The Christoffel Darboux theorems of section l^?^ show that the kernels can be computed 
in terms of only with n — d2 < m < n and with n — l<m<n + d2 — 1 (resp. 

with n — di < m < n and ^/'^ with n — 1 < m < n + di — 1). We thus introduce 
the following vectors, called "windows" P ITIH ITT]: 

ll)n{x) = {iJn-d2,i'n-d2+l, ■ ■ • ,^n-l, V^n)* = n""'^2V'oo (6-1) 

and similarly: 

4(2/) = nr^^^^^ , V^Jy) = n::;^^_,^"oo , ^Uy) = K-X_,<f^ (6-2) 

We also introduce the following matrices, called "windows" HUl E] : 

'^nix)mi ■= , m = n - (^2, • • • , n , i = 0, ...,^2 (6-3) 
'^'n(2/)mi := 0^^^!/) ) m = n-di,...,n , i = 0,...,di (6-4) 
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^n{y)mi ■= "tPmiy) ) ITl = U - 1 , . . . , U + di - 1 , i = 0,...,di (6-5) 

^n{x)mi ■■= (f>^r^{x) , m = n - 1 , . . . , « + - 1 , i = 0,...,d2 (6-6) 
The relationship Eq. ()4-H()|l imphes the duahty (found in [Oj): 

$J,(x)A„*„(x) = Id (6-7) 

and similarly: 

<^MB'Jn{y) = ld,,+, (6-8) 
Notice that (combining Eq. ()3-30|) and Eq. ()6-7|) ): 

^oo(x) = (L(x) - i?(x))A„^„(x) (6-9) 

The matrix 

Fn{x) := {L{x) - R{x))An (6-10) 

is the so called Folding matrix of ^T]. It's property is to fold any operator acting on 
\E'oo into an operator acting on the window only. For any finite band operator O, we 
have: 

K'"' o^oo{x) = (nr^^ 6i^„(x)) ^„(x) (6-11) 

n^~'^2 OF„(x) is a square matrix of size + 1 x ^2 + 1. When acting on the window 
\E'„, it gives the same result as the operator O acting on the infinite matrix \E'oo. 

7 Mixed correlation function 

We now arrive to the main results of this article, which concerns the computation of 
the mixed correlation function 

Wn{x, 2/) = 1 + / Tr \, \, ) , (7-1) 
\ X - Miy - M2/ 

where Mi and M2 are in the ensemble iJ„ x HniV) of remark [2.11 and explained in 
[TS] . with the measure: 

and where Z is the normalization constant, called partition function: 

Z=f e-^iViiM^)+V,iM,)+M,M,) ^7_3) 

We recall that iJ„ x i/„(r) is the set of normal matrices (i.e. [Mi, M|] = = [M2, Mj]) 
with pairs of eigenvalues constrained to be on F, see ^H] for more details. In the case 
F = R X R , Ml and M2 are hermitean matrices of size n. 
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The formula of [H] reads: 



(7-4) 



det Id„ + n„_i 



x-Qy~P' 



■n„_i 



which is a n X n determinant, and is thus not convenient for large n computations, and 
for many other applications. 

The purpose of this section is to write it in terms of a determinant of size d2 + l or 
di + 1. 



7.1 Formula for Wn 

We introduce the following lower triangular matrices^ of size di + 1 (resp. d2 + !)'■ 



Un{x,y) 
Unix,y) 



y + vi{x) 

In 



p p* — TT 



" + ^^'^^Ue* -n 



X - Q 



(7-5) 



7n 



y-P 



Many of their properties are described in appendix C. 
The following theorem is proved in appendix D: 



Theorem 7.1 



Wn{x,y) 

-flKn+i{x,y) Jn-i{x,y) det(Id„ - n„(l - 



n+l 



Jn-1 



7lK^+,{x,y) Jn-iix,y) det(Id,,+i - U^-'^il - |^)f/„(l - ^W^] 



n+l 
t 



Jn-1 



llK^^,{x,y) Jn-i{x,y) det(Id,,+i - U^-'^^il - ^p^l - ^)U^) 



(7-6) 
or also: 



Jn-1 



Kn+l ' 



Wn{x,y) = -f'^ det {Id n+l -Unix, y)Un{x,yy) 

{Jn+d2 + 4lZ^^n) {Kn-d2 + 0n i_J„^t V^n) 



(7-7) 



They are infinite matrices, with only a non- vanishing sub-block of size di + I (resp. d2 + I)- 
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Notice that those formulae involve only functions which are within the windows. The 
determinants are in fact of size mm{d2 + l,di + 1). 

The proof is given in appendix D. 
Theorem 7.2 

Wnix,y) = -i Kn^d2ix,y) Jn+d2ix,y) det{y - M{x,y)) (7-8) 
where M is the (^2 + 1) x (c?2 + 1) matrix 

M., = Hl:^\x,x) + (y - d^)K^:il{x,y) + '^^ (y + d^)J^Z^{x,y) 

(7-9) 

proof: 

We use expressions of appendix F, and Lemma Ol □ 
7.2 Recursion Wn+i — Wn 

We introduce the following lower triangular matrices of size di (resp. ^2): 

>v„w^n„K(£)^n» . >v„fe) ^ n„liM^n" (7-10) 

X - Q y - 

The following theorem is proved in appendix E: 
Theorem 7.3 

Wn+i{x,y) - Wn{x,y) 

t ' * 



(7-11) 



J„ det(l - n„(l - %^)W„(1 - ^^)>V*) 



I.e. 



Wn+i{x,y) - Wn{x,y) 

det(l-W„W*) . [[J^ + {^^Jwij^^Wj^mJ^^^^^) 



(7-12) 

i.e. Wn+i — Wn actualy involves the computation of a determinant and inverse of a 
matrix of size min((ii, ^2). 

The proof is given in appendix E. 
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8 Application: differential systems and spectral 
curve 



An important observation of [0] and pTH |321 ISH] is that windows of consecutive bi- 
orthogonal polynomials satisfy some integrable differential systems, and a Riemann- 
Hilbert problem. An explicit representation of those differential systems was found in 
[TT] . Here, thanks to the result for the mixed correlation function, we are able to give 
new representations of those systems, and compute explicitely their spectral curve. 



8.1 Differential systems 

We define the following (^2 + 1) x [d2 + 1) and {di + 1) x [di + 1) matrices: 

VM := Kiy) K\y) , ^n(y) := K\y)- 



(8-1) 
(8-2) 



Those matrices are square matrices of the size of the corresponding window, and they 
have polynomial entries. They give some ODE's for the windows: 



<(x) =r'„(x) ^„(x). 



-3) 



It was found in [9 that they enjoy some duality relations (which merely come from 
duality Eq. (jEH)): 



Vl^{x)An = AnVnix) , V^{y)Bn = BnVn{y). 
Theorem 8.1 We have the following equivalent expressions for T): 



-4) 



Vn{x) 


= Ill-'^^P\L{x 


) - R{x))An 


y - Vn{x) 


= {Ui-\x,y)- 


Ul-^^Un{x,y))An 


(vl/-i(x)l?„(x)vl/„(x)),^. 


= H^'^\x,x) 





(8-5) 

And the matrix elements of y — Vnix) are {di -|- 1) x [di + 1) determinants given 
by (0<k,l< d2): 
-Ifk>di 



det 



{y - Vr,{x))^_^^^_i = 

1 lm=l ^n— 171,71— m—d2 



X 



nrf {x - Q)K--Z-T nrf {x - g)e„_, 
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(8-6) 

-Ifk<di 

(y - ^^n(x))„_,,„_, = ^~i.r ' X det 

llm=0 In+m—l [lm=l ^n-m,n-m—d2 

(8-7) 

The first equality of Eq. ()8-5|) was found in ^1] , it comes from the Folding matrix 

Eq. (i?mn) : = n;^-'^^^'^ = xii-'^^p'^^ = u'^-'^^p'f^{x)'^^. 

The second equality is a rewriting of the first one using the matrices Un and Un, 
see appendix C. 

The third equality is a mere rewriting of the definition of V, using the duality 
Eq. (0301). 

The last representation is obtained by using the first formula, and computing the 
folding matrix by inverting the linear problem \l/oo = -^n^n in the kernel of x — 

8.2 Spectral curve 

It was proven in Pj that all those differential systems share the same spectral curve. 
We define: 

£nix,y) := t det(?/-T'„(a;)) =t det(?/-P„(a;)) 
tdet{x-Vniy)) = tdet{x-'Dniy)) 

(8-8) 

We are going to give here several equivalent formulae for computing S. 
Theorem 8.2 We have: 



£:„(x,y) = -7^det(l-f/„?7*: 

(8-9) 
and 

£n+iix,y) - £nix,y) = -det(l - W„W*) (8-10) 

proof: 

For the first expression, use the 2nd expression in theorem 18.11 The proof of the 
recursion formula is found in appendix H. □ 

Then, we introduce the following Lemma, which consists in taking the polynomial 
part at large y of formula Eq. ()7-7|) : 
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Lemma 8.1 



(8-11) 

^^^^ \ x-Mi y-M2 I iW ' ^ ^ 

Notice that this formula is the finite n counterpart of what was found in the formal 

large n expansion in ^1], namely theorem 3.1, equation (3.3) of [Tlj . 

proof: 

This Lemma is proved in appendix F.5. □ 

Using this Lemma, taking the polynomial part in x at large x, we prove the following 
theorem, which was conjectured by Marco Bertola [Sj: 

Theorem 8.3 Proof of Bertola's conjecture: 

Pol Pol {V;{x) + y)iV^{y) + x)Wn{x, y)=2n + £n{x, y) (8-13) 

X^OO J/— >oo 

i.e. 



(8-14) 

Marco Bertola has proved it for potentials of degree ma.x{di,d2) < 5, which is quite 
remarkable, and also for the smallest values of n 

This theorem is also the finite n counterpart of theorem 3.1 of ^3]. Notice that 
the spectral curve for finite n is the same as the algebraic curve found from large n 
considerations This fact has important consequences, some of which are described 
below in the next section, 
proof: 

The proof is in appendix G. □ 

For completeness, we also write a formula for the spectral curve, due to Jacques 
Hurtubise: 

Theorem 8.4 Hurtubise formula \2(^ 

Six, y) = ^ det ( " y\ \ (8-15) 
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which simply amounts to say that if Sn{x,y) = 0, there must exist some functions 
ipn-d2~i: ■ ■ ■ yi^n+di, such that we have simultaneously yipi{x) = ip'i{x) = J^i^jii^ji^) 
for i = n — d2, . . . ,n and xtpilx) = Qijipj{x) for i = n — l,...,n + di — 1, i.e. the 
vector {%l)n~d2~ii ■ ■ ■ 1 "^n+diY IS the kernel of the matrix above. 

9 Examples 

9.1 Gaussian case, Ginibre polynomials 

Consider Vi{x) = V2{y) = 0, which is Ginibre's ensemble. It is well known that the 
bi-orthogonal polynomials are monomials 



Pn{x) = , qn{y) =1/" , hn = nlvT, (9-1) 
and the Cauchy transforms are 

vrn! Tcnl 

P-(y) = ^ ' ?n(x) = ^. (9-2) 

We have: 



Qnm = Pnm = Vn + 1 5m,n+l , 7n = V^- (9-3) 

We find that for n > m: 

Rnm{x) = -V^„(x)0m(x), (9-4) 

and is an upper triangular matrix. 
In that case, theorem 17.11 becomes: 

Wn{x,y) = l + nJn-i{x,y)Kn{x,y) - xyJn{x,y)Kn-i{x,y) 

= nKn+iJn-i - xyKnJn, (9-5) 

and theorem ()7.3p with W„ = becomes: 

Wn+l -Wn = K^+lJn. (9-6) 

9.2 Gaussian elliptical case 

Consider Vl{x) = tx, ^'(i/) = ty, we write: 

6 = l-ti (9-7) 
The bi-orthogonal polynomials are rescaled Hermite polynomials , and we have: 



O'ri Qn—l,n Pn—l,n \l ^ (9 
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Qn,n—1 t^n 

The Christoffel-Darboux matrices are: 

1 



An = 7n 

The matrices f/„, f/„ are: 
Un{x,y) = - 



t 



n,n — l 



Bn = 7ri 



1 
t 



t 

7n 



In 

i 



and the matrices VV„ and Wr, are of dimension 1: 



W(a;) = t 



The differential systems are: 



thus the spectral curve is: 



t t 

'57n -ty 



VJx] 



i>n{y) 



-tX S'Jn 

t i 

-iy 6jn 

t in t 



£n{x, y) = idet{y - P„(x)) = (x + iy){y + tx) - 6n = -7^ det(l - U^U^J 
Theorem 17.31 gives: 

W^n+l ~ = JnKn+l — ttJn+lKn 



10 Tau- function 

Theorem 18 . 31 implies for k < di and I < (12- 

{ Tr Ml) = Res Res ^"i^'wl?.'^' rfy 
Thus, if we introduce the matrix integral: 

Z= /'dMidM2e-EM*MTrAf^4] 



we have: 



S^iix ^ y^ x^ y^ 
Res Res ", , \ -, . dx dy 

x^oo y^oo V{{Xj ^2(2/) 



dlnZ 



dtk. 



tk.i=0 if /£/>! 
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Equation Eq. ()10-3p is similar to the tau-function of [2HI HZl IHH 122] generalized to the 
resonant case by 7 , and thus, similarly to 12], we have proved that: 



T{{tk,i}) = Z hih,i = Oifkl>l. (10-4) 

This strongly suggests that the general setting of the isomonodromic problem should 
include all tk/s. 

It is easy to define bi-orthogonal polynomials with a weight of the form 

uj{x,y) = e-^k,iik,ix''y' 

and they are related to the so called ensemble of normal complex matrix integrals (see 

^normal = / ^Mc-^m^m-m^m. 

where Nn is the set of complex matrices which commute with their adjoint [M, M^] = 
(they have complex eigenvalues not constrained on any path). For normal matrices, 
we have: 

Tr — i ^\ = Tr ( n„_i — ^ , ^ n„_i | (10-6) 

so that our Wn is not the mixed correlation function for normal matrices. 
Our Wn{x,y) is in that model: 

V) = (det (l + ) = det(l + n„..^^n„_0. (10-7) 

It is not clear yet why this function should be related to the spectral curve for the 
normal matrix model... 



11 Conclusion 

In this article, we have found a formula for computing the mixed correlation function as 
a min((ii + I,d2 + 1) determinant, instead of a n x n determinant. This kind of formula 
can be very useful for finding large n limits. Beside, we have found several expressions 
for the spectral curve, which can be useful for studying the integrability properties. 
In particular, we have proved the conjecture of Marco Bertola, which should open the 
route to other applications. 

In this article we have computed the 2-point mixed correlation function, however, 
for many applications to physics (BMN limit of ADS/CFT correspondance in string 
theory, or boundary conformal field theory), one needs mixed correlation functions 
involving more than one trace, and more than two matrices in a trace. The key element 
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of the work of [S], was the use of Morozov's formula for 2-point correlation functions 
of unitary integrals. Since then, some generalizations of Morozov's formula have been 
found for arbitrary correlation functions [21], and it seems natural to generalize the 
result of [HI, and of the present paper. In that prospect, one could expect to understand 
the Bethe-ansatz-like structure of large n mixed correlation, which was found in [20] ■ 
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Appendix A Proof of theorem 13.11 

proof: 

Notice that ipniv) is piecewise analytical in the connected domains separated by 
the ^/^'-^'s. In each such domain, it bahaves at oo, like: 



My) ~oo V^r"-' e^^(^) (1 + 0{l/y)) (1-1) 
Let fn{x, y) be the primitive of ipniy)^^^ which vanishes at oo in sectors where KeV2 < 0: 

fn{x,y)= F My"V'" (1-2) 



fn{x,y) behaves at oo in each connected domain, like: 

^ oo 

fn{x, y) -oo y vd/-"~'^~' e^^(^)+^^ (1 + ^ /„,,(x) y') (1-3) 



k=l 



where each fn,k{x) is a polynomial in x. 
By definition of Lnm{x) we have: 



TT-T] dy<Pm{y)e {fn{x,y+) - fn{x,y-)) 

2«Vr ^ y;r,(,) 

77— y] / dy(p^{y)e-''y fn{x,y) 

Im Jar) 
domains D 

domains D "^^-^ fc=l 

00 

= Res dy qUy) y'^-'^-' e"^^ (1 + V /„,fc(x) y-'^) 

fc=i 

(1-4) 

This expression vanishes if m < n + (^2, and is clearly a polynomial in x. 
Now, compute: 



d2 



e-v'2(y)gV2(y")e-^2(?/')e-^i(^') 



1=1 

] V / dy [ dy" I dx'dy' 

2i7ry/hnhm "^Aw if'' ir 

] 7, } Pn{x )qm{y) e ^ e ^e ^ 

X — x' y" — y' 



m J"f(^'> Jf^'' J^Cfe) J^O) 
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1 Vi{y")-Vi{y') 



X — x' y" — y' 



d2 



= EE 



2m \/hnhr, 



•"m J 7 



[ dy i 



dy" / dx^ / dy' 



(k) 



i=l j,k 



(1-5) 
Moreover we have: 



X — x' y" — y' 



— f dy" [ dy' e^-iy")e-^-^y') ^^^y"^ ~ ^^^y'^ e^y" e^'y' 

— / dy" I dy'Y^. 

217: Ao) y" - y' 

2?^ .^«(«iu^Wf«) Ao) y"-y' 



y - y 



- / dy" f dy'- — 

™ J¥+{yW-{y) A'^' yy" - y' 

- I dy" [ dy'—^, 

^'^ J¥i\v)u^^\v) Ao) y -y 



.QV2{y")^-V2iy') ^-x'y' ^xy" 



2m 

2m 

1 



eV2(y")g-V2(2/') ^-x'y' ^xy" 



dy' 



yll _ yl 



.QV2{y")Q-V2{y') g-x'y' ^xy" 



l''(3/) [y -y) 



.QV2{y")Q-V2{y') Q-x'y' ^xy" 



dy" I dy' _^e^^(^")e-^^(^') e"^'^' e^^" 



= 5ij Res rf?/' ■ 



y-y 



-- / 

2m M 



'r+'iy)ut:'{y) J J 



(1-6) 



Therefore: 



d2 



2^7r ./^«(2;)uf«(2/) ^7 



y"- 


-y' 


-x'y' 


^xy 


X — 


x' 


y"- 




-y' 


X — 


x' 



- [ dy" [ dy' ^^e''^^y"h-^'^y'^ e-^'y' e^y" 

2?7r i-W(„^u^«f„) 7^0-) y"-y' 



y - y 



d2 



= -EE 



i=l j,k 

d2 



-EE 



Hkj 



yd) 



dy / dx' e-^^(^)e-^^(^') -Pn{x')qm{y) e"^' 

J-y(k) X-X' 



=1 j,k 



2m \/hnh. 



rn J 7' 



dy 



dy" / dx' / dy' 



(fe) 



0) 
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y - y 

1 



X — O , 



yy 



- Q ; 2z7r J-w (j^)ufj) (j^) 



x-Q 

Using Eq. ()2-52p . this also proves that L is a left inverse oi x — Q. 
□ 

Appendix B Proof of determinant al formulae 

B.l Inverses 
1. 



n„-i 1 - n„(i - ^)n„_i 

"^'^ V V'n V'O V'n V'O / 



n„_Mi-^in„_. 

Id 



n 



2. inverse of n„_i^n„_i: 



^-Q J \ (f)n-l 

iin-i{R + Voo0oo n„-i)(i + n")(x - Q)n„_i 

{x - Q)Tin-l 
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= n„_i(i? + ^oo0oo*(n„_i + n"))(x - g)n„_i 

= Idn (2-2) 
3. inverse of n„_i(x - Q){y - P*)n„_i: 

(n„_i(x -Q){y- P*)n„_i). (n„_i^*n„ (^i - -^Vooc/'oo*) n„i?n„_i) 
n„_i(x - Q)(y - p*)i?*n„ (i - n„i?n„_i 
n„_i(x - g)(n„ - ^n„) (i - n„pn„_i 
= n„_i(x - g)(n„ - - + ^n„^)n„i?n„_i 

n„_i(x - Q)(n„ - - + ^)n„it:n„_i 

n„_i(x - g)(n„ - n„^^n„)i?n„_i 
n„_i(x - g)(i - ^^n„)i?n„_i 
n„_i(x-g)i?n„_i 

Id„ (2-3) 



4. • Proof that: 



(2-4) 

We have: 



a; - Q 2/ - P* 

rin-l 7Tn„_2 — -n„_i -l- Jn-l^n-l4'oo4'oo^^n-\ 

x-Q y - 



1 



nm 



X - Qy - P 

If 11 
/j—r— / Imiy') ; -pn{x') uj{x\y')dx' dy' 

y/ tin rim Jr y ^ y ^ 

I cLiJ,[x , y 



VKhm Jr y-y' X- x' 

rr—r— / -f uj{x ,y)dx dy 

vhn hm Jv y - y X - X' 

/ ) 7 u;{x,y)dx dy 

Jr y-y x-x' 



Vhn h. 



m Jr 

1 f Qm{y) Pn{x)-Pn{x') , ^ f ^ > 

uj{x ,y ) ax dy 



VKhm Jry-y' x-x' 

1 / Qmiy) - Qmiy') Pn{x) - Pn{.x') , ,s , / , , 

+ nr-f- / zr—:, — ^{x,y)dx dy 

yii'nii'm Jr y y x — x 
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i/jn{x)(l)m{y) Mx,y) 



vhn hm Jry- y 

n—l m—l ^ „ 

+ V V rr-^ / Rmiiy) Rnjix) qi{y')pj{x') uj{x\ y') dx' dy' 

= (^Jo{x, y) ipoo{x)(t)oJ'{y) + 'ipoo{x)(t>J {x)R\y)+ 

R{x)^^oo{y)<pJ{y) + R{x)R\y)) (2-5) 

/ nm 

Thus: 

^— ^ = ^0 i^oo(t>J + ipoo(fjR' + Rifoo(f>J + RR' (2-6) 
Now, we take the n„_i projection: 

= n„_i'0oo0oo*-R*n„_i + n„_ii?'0oo0oo*n„_i + n„_ii?^*n„_i 
= n„_i'0oo0oo n„_ii?*n„_i + n„_ii?n„_i'0oo0oo*n„_i+ 
n„_ii?n„_ii?*n„_i 

= n„_i(i? + -0000^00 )nn-l(^* + '0^oo0oo*)nn-l 

-n„_iV'oo0oo nn_iV'oo0oo n„_i 

= ( Jn - Jo) n„_iV'oo0oo*n„_i 

+n„_i(i? + '0oo0oo*)nn-l(^* + '0oo0oo*)n„_i 

= n„_i — T^nn-i — — n„_i + j„n„_iV'oo0oo*n„_i (2-7) 
X - Q y- P* 

Similarly we have: 

rin-i^ — — ^n„_i — Jo n„_i^/'oo0oo n„_i 

= n„_i'0oo0^oo*^*n„_i + n„_ii?^'^oo0oo*n„-i + n„_ii?^*n„_i 

= n„_iV'oo0oo n„_2^*n„_i + n„_ii?n„_2^oo0oo*n„_i + n„_ii?n„_2^*n„_i 

= n„„i(i? + v^oo0oo )n„_2(-R* + v^oo0cx)*)n„_i - n^-iV^oo^oo n„_2^oo0oo*n„_i 

= {Jn-l - Jq) n„_i?/;oo0oo*n„_i + Tin-l{R + V'oo^oo )n„_2(-R* + i)oo(t><J)^n-l 

= n„_i — 7:n„_2 — — n„_i + j„_in„_iV^oo0oo*n„_i (2-8) 

x-Q y - P* 

5. • inverse of n„_i^^^n„_i: 
First, notice that: 

(fjir'-\x - g)n„-i = (fj{i - ^n-2){x - g)n„-i 
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= 4>oo {x - (5)n„-i - n„_2(x - (5)n„_i 

- (/);o*n„_2(x - Q) 

Wo 

(2-9) 

Now, compute: 

Jn-\ ^ 

n,_i(7/ - P*)(n,_2(a; -Q) + -^n"-V;o(T^ - 0OO n„_2(a; - Q))) 

= n„_i(y - P0((i - n"-^^^)n._2(a; - g) + n'^-ii:2^) 

(2 - 10) 
Thus: 

- t 

n„_i(y - p*)(n„_2 + n"-i^^^n"-i)(x - g)n„_i^^— n„_i 

Jn-i x-Qy - 

= n„_i(y-p*)(i-n'^-^^^)n„_2^-n„_i 
+n._i(y - P^n-^-^e^ \ ^ n._, 

J„_iV'o x-Qy-P^ 

= n„_i(i/ - p*)(i - n--i^^;^)n„_2(p* + i^^<pj)iin-i 

+n„_i(y - P*)n"-l-^e* (P + V'oo0oo>n-2(^* + V'^oo0oo*)n„_i 

+n„_i(y - p*)n«-i-^e* (j„_iV'oo0oo')nn-i 

Jn-lWO 

= n„_i(y - pO(i - n«-^^^^)(p* + iin-2^foo^J)Tin-i 

+n„„i(|/ - P*)n"-1^0;o*(i?* + Iln-2i^oo(t>J)^n-l 
Jn-l 

+n,_i(i/ - p*)n^^-v;o0oo*n„-i 
= n„_i(y - p*)(p* + n,_2^;o0oo*)nn_i + n„_i(y - p*)n"-V^oo0oo*n„_i 
= n„_i(y-p*)(P* + V''oo0oo*)n„-i 

= Id„ (2-11) 
.2 Determinants 

1. • determinant of n„_i(a; — g)n„_i. This is a classical result and can be found 
in the litterature. 
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One possible proof, is that det(n„_i(x — (5)n„_i) is a monic polynomial of degree 
n. Consider x a zeroe of this implies that Iln-iipoo{x) — Iln'ipoo{x), thus: 

n„_i(a; - Q)nn-iipoo{x) = n„_i(x - Q)nnip^{x) = n„_i(a; - Q)tlJoo{x) = 

(2-12) 

Thus, the n zeroes of V'n are also the zeroes of det(n„_i(x — Q)Iln-i)- The 
converse is easy too. 

Thus we have: 

det(n„_i(a; - g)n„_i) = pn{x) , det(n„_i(y - P)n„_i) = (2-13) 

2. • determinant of Un-i^zQlin-i' 
Compute: 

p„(x) det(n„_i — 3— Iln-i) 
a; - 

= det(n„_i(x - (5)n„_i — 7:n„_i) 

X — Q 

= det(n„_i(x - g)n„_i(i? + v^oo^oo )n„_i) 

= det(n„_i(a; - Q)(n„ - e„e^)(i? Voodoo )nn-i) 
= det(n„_i(a; -Q){R + V'oo^oo )nn-i 

-n„_i(x - g)e„e^(i? + -000000 )nn-i) 

= det(Id„ 7„en_ieJ,(i? Voodoo )nn-l) 

= 1 7„e^(i? + V'oo^oo )e„_i 

= 1 + 7n(-Rn,n-l + V'n^n-l) 
= l + 7n( KV'n0n-l) 

= TnV'n^n-l (2-14) 

V V 

det(n„_i^-n„_i) = ^^^^ = ^== (2-15) 

3. • kernel K^- 

det(n„_i(a;-Q)(|/-P*)n,_i) 
= det(n„_i(a; - Q)U^-i{y - P*)n„_i + n^_,{x - Q)U^{y - P')U^-i) 
= det(n„_i(a; - Q)n„_i(y - P*)n„_i + TX-ieJ^-i) 
= Pn{x)qn{y) det(Id„-F 

+7^n„-i f 1 - n„Pe„_ie* _,p*n„ f i - n„_o 



Pn{x)qn{y) det(Id„ + n„_i (^1 - ] e,,el ( 1 - ] n,_i) 
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P.M(l + <^(l-^)n„.,(l-*|<)e„) 

Pn{x)qn{y) (1 + — ^) 
Pn{x)qn{y) 

h^ey^i.^)+y^^y) Kr,+i (2-16) 



4. • kernel J„: 



det(Jnn„_i^/'oo0oo*n„-i + n„_i — 7Tn„-i — ^^n^-O 

a; - Q y - -P 

(1 + j„0 

oo (nrt— 1 

2; - V y - -P 
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(1 + J„0oo*(n„_i^— n„_i)-i(n„_i^— n„_i)-Voo) 

y - x-Q 

g-Vi-Va ^ 

— Vn-10n-l(l + Jn0oo*n„_i(?/ - P*) 

(1 + n"^^:^)n„_i(i + ^^^^n")(x - g)n„_iV'oo) 

— V'n-10n-l(l + ^n0oo*n„_i(|/ - P*) 

(1 + n"^i^^)n„_i(i + ^^^^n")(a; - g)n„_iV^^) 

^n-l 0n-l 

g-yi-V2 ^ 

— ?/'„_i0„_i(l + J„0oo*n„_i(y - P*) 

t 

(1 + n"%^)(n„_2 + e„-ie*„_0(i + ^^^=i^n-)(x - Q)n„_iVoo) 

^n-10n-l(l + Jn0oo*n„_i(y - P*) 



(1 + n":^^^^)e„_ie:,_i(i + ^^^^^^ii-){x - g)n„_iV'oo) 

Wn-l 4>n-l 

— V;n_i0„_i(l + J„0oo*n„_i(|/ - P*) 

(e„_i + n"i^)(et,_i + |^n")(x - g)n,-i^oo) 

tPn-l (Pn-l 

Q-V1-V2 ^ 

— V'n-10n-l(l+ 
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-V2 


hn- 


-1 




-V2 


hn- 


-1 




-V2 


hn- 


-1 




-V2 


hn- 


-1 



fn-l (f>n-l 

- t 

1pn-l (Pn-1 

'0„_10„_l(l + Jn(t>J^ —-^ ^^oo) 

'0„_10„_l(l + ) 

Jn-1 (2-17) 



Appendix C Matrices U and U 
C.3 Definitions 

The following matrices play an important role. 
Define: 

Unix, y) := Univ - P')Rix)U^-^ , Unix, y) := n„(x - Q')Riy)U"-' (3-1) 

Notice that since R is lower triangular and P is finite band, Unix,y) (resp. Unix,y)) 
is a small lower triangular matrix of size di + 1 (resp. d2+ 1): 

Unix, y) = nr '^^ iy - P')Rix)K+i-i (3-2) 

Unix, y) (resp. Unix, y)) is linear in y (resp. x) and of degree at most di in x (resp. ^2 
in y). 

Unix,y) can also be rewritten as: 

U^i., y) = -l^±llMe.e:.,, - n„liM^n"-. (3-3) 

proof: 



c/„(x,y) = nr'Hy - ^*)^(^)n;:;i_i 

= nr'^Hz/ + K(^) - viix) + y/(g) - y/(g) - p')Rix)Tii-\_, 

= nr'^^y + Viix))Rix)Ill-\_, - KT'^Viix) - ViiQ))Rix)Ill-\^_, 

-K-''iviiQ) + P*)Rix)Kl\-. 
= K-''iy + ^i'(^))^(^)n::;i_i - K-Hviix) - viiQ))Rix)xiiA-. 

= iy + Viix))Ill-''^R{x)Tll-\^_, 
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7„ e„e;,_, - \;_g (3-4) 
□ 

Notice also that if n > (^2^1, we have: 

Ur,{x,v) := IlnR{x){y - P*)n"-' (3-5) 

C.4 Some properties 

• Multiphcation by the Christoffel-Darboux matrix: 

- Un{x, y)Ar, = + {y- P*)ene^ - n„(y - P*)Pn"(a: - Q) (3-6) 

- A^C/^a;, I/) = P*_i + e,_iej,_i(|/ - P*) - (x - Q)n„_i(|/ - P*)Pn"-i (3-7) 
and subsequentely: 

¥jJn{x, y)Ar,^ 00 = -C^n+l^oo (3-8) 

^'^Ar,Un{x, y)^oo = -l>L5^-l*oo (3-9) 

proof: 



-Un{x,y)An 

n„(y - P*)P(a; - g)n" - n,(|/ - P*)Pn"(a; - Q) 

n„(y - p*)n" - n„(y - p*)it:n"(x - Q) 

[n„, y - p*]n" + (y - p*)e„e^ - n„(y - p')mr{x - Q) 

= Bi^, + (y - P*)eX - n„(y - P*)Pn"(x - Q) 
(3-10) 



-AnUn{x,y) 

iin-i{x -Q){y- P*)pn"-^ - (x - Q)iin-i{y - P*)i?n"^^ 
n„_i(i -f (y - p'){x - g))pn"-i - (x - g)n„_i(t/ - p*)pn^ 
n„_i(y - p*)n"-i + lin-iBW-^ -{x- Q)Tia-i{y - P*)pn" 
n„_i(t/ - p*)n"-^ - (x - g)n„_i(y - p*)pn"-i 
n„_i[i/ - p\ n"-i] + n„_in"-i(2/ - p*) 
-(x-g)n„_i(i/-p*)pn"-^ 
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Iin-i[P\ n„_2] + e„_ie^_i(|/ - P') - (x - Q)Iin-i{y - P')Bir-^ 



(3-11) 

□ 

Inverse: 
proof: 



ui-\x,y) = K-'\y - P')mKll-. (3-12) 



t>n(^,i/)nr'^^(i/-P*)i.(x)n::;^^_i 
= n-ip*(i/)(a; - g)nr'^^(i/ - P')L{x)Iil-\^_, 
= n"-ip*(y)n-(x - g)n,(?/ - p*)n„_iL(x)n:j;^^_, 
= n"-^p*(y)n"(x - Q){y - P'nn-M^)Kl\,^. 
= n"-^p*(y)n"(y - p*)(x - g)n„_iL(x)n::;^^_, 

= Tr-'R\y){y - P'){x - Q)L{x)Xll-\^_, 

= K-+\,-i (3-13) 

□ 

Multiplication of the inverse by the Christoffel-Darboux matrix: 

-Ui-\x,y)A^ = Bi_,jr + {y-P')er,el-ni-''-{y-P')L{x)Tr{x-Q) (3-14) 

AjJl-\x,y) = -W^Bi^,^ + K+a._,{y - P*) - (x - g)n«(|/ - P%{x)Yil-l^_, 

(3-15) 



In particular: 



and 



^'^ij'-\x, Z/)^„*oo = -$L^n-d,nn-l*oo (3-16) 

^'^AjJl-\x, y)^ = -|.*^n"P^^.,^^ (3-17) 



proof: 



Ul\x,y)A^ 

nr'^^(y-p*)L(x)[a:-g,n„_i] 

nr'^^d/ - P')L{x){x - Q)Iin-i - K~''{y - P')L{x)Y[r.-i{x - Q) 

Tll-'^^{y - P*)n„_i + -fn-d2Ln-d2-l,n-ien-d^^i-l{x - Q) 
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(3 - 18) 



= K-'Ky - P')L{x){x - - K-'^iy - P'WxWix - q) 

[K-'^\y - p*]n" + {y- p*)n;^-'^2n" 
-nr^^(y-p*)L(x)n"(x-g) 
= p*_,^n" + (y - p*)e„e:, - nr'^Hy - p*)L(x)n"(x - g) 

(3-19) 



Then: 



AnUi'\x,y) 

-[x-Q,W']iy-P%ix)Ul-^ 

n"(x - g)(i/ - p')L{x)ui^l_, - (x - g)n'^(i/ - p*)L(x)n::;i^_, 

U-il + iy-P')ix~Q))Lix)n--l^^, 



jn-l 

-d2-l 



-{x-Q)U-{y-P%{x)K-X_, 
my-P')K-i_, + U-Lix)K-X_, 
-{x-Q)U-{y-P%{x)U--\^_, 



u-{y - p*)n„+,,_, - (x - Q)u-{y - p*)L(x)n;:;i^_, 
n"[(y - p*), n„+,,_i] + n::+,^_i(y - p*) 

-(x-g)n"(y-p*)L(x)n::;i^_, 
-n-p*^,^ + n::+,^_,(y - p*) - (x - g)n'^(y - p^ix)^-^^. 



(3 - 20) 

□ 



Appendix D Proof of theorem 17.11 

proof: 

Let us introduce the following n x n matrices: 



Cn:=n„-i(x-g)(y-p*)n,-i (4-i) 

:= nn-i^^n„_, (4-2) 

Un := 1„ + n„_i(|/ - P')Un-2{x - Q)U.n^l (4-3) 
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and, according to formulae of section 15.11 

detCn = /i„i^n+i , c' = n„-ii?*n„ (^1 - -^tAoo^oo*) n„i?n„_i (4-4) 

deten = ^ (4-5) 

= n„_i(y - p*)(n„_2 + -^n"-V"oo0oo*n"-i)(x - Q)n„_i (4-6) 

We have: 

= i„ + n„_i(?/-p*)n„_2(x-Q)n„„i 

i„ + n„_i(y - p*)(x - Q)n„_i - n„_i(y - p*)n"-i(a: - g)n„_i 
= Cn - n„-i(y - p')u"-\x - g)n„_i 



1 



(4-7) 
Thus: 



C„(l - -n„_i(y - P')U^"\x - Q)U 



n-1 



(4-8) 

and: 



det a;„ 

det C„ det(l - ^U^-i{y - P')U^-\x - Q)n„_i) 

KKn+1 det(i - n„_iP*n„ - -^^v^oo0oo*^ n„Px 
x(i/-p*)n"-i(x-g)n„_i) 

det(i - n„_iP*n„ - -^^^Aoo^oo*^ f/n(a; - (5)n„„i) 

/i„ir„+i det(i - n„ ( 1 - -^^^<pj\ Uniil) 



detUn = hnKn+i det{l - UnU^ + — — U^^^cpJUnU^) 

= KK^+, det(l - f/„t/*) (1 + )^-> c 
K det(l - f/,t/*) {<PJ- ]-j-i^oo) 

(4-9) 



Inverse: 
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(1 - ;^n„„i(y - p')u-'\x - Q)n„_i)'' C 

(1 - n„_ii?*n„ (i - ^^1^) u^ix - g)n„-i)-' 

(1 + n„_ii?*n. ^ — - — ^ — (i - ^^"j ?7„x 

1 - n„ (i - j f/„f/* V ^n+i ; 

(x-Q)n„_i)c' 

+ n„-ii?*n„ — — ^n„, (i - 1^^] f/„x 

1 - n„ ( 1 - f/„[/* V ^n+i ; 



n„_i^*n„ 7 — ^——^ — - fi - I n„i?n„_i 



1 _ n 1 - 1^^^ u ijt V -^n+i 



(4 - 10) 

Now, the formula of [H] gives Wn = det(l + ^ri), thus we compute: 

det(l + a) 

det^n det(l + C') 

deta det(i + n„_i(y - p*)(n„-2 + -^n"-V"oo0oo n"-i)(x - g)n„_i 
det^n detK + -^n„_i(y - p*)n"-v;o0oc n"-i(a; - g)n„_i) 

deta detcu„ (1 + -^cfjw\x - Q)Ii^^^uj~^Iin_^{y - P*)n"-V"oo) 

(4-11) 
Thus: 



n 

t 



detu;„ (^-1 + 0ooV* J ^ s :- ( 1 - ] f/„^„ 



(4 - 12) 



1 _ n 1 - r/ r/* v -'^"+1 



I.e. 

W^n = 7^i^n+i^-i det (1 - n„(l - %^)f/„(l - ^^^)f/* I (4-13) 

The two other formulae are obtained by Weinstein-Aronstein duahty. 
The lemma Ol gives Eq. ()7-7p . 
□ 
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Appendix E Proof of theorem 17.31 

The formula of jH] gives: 



Wn{x,y) = det Idn + TTn-i 



n 



1 1 



x-Qy-P 



J 

n— 1 



(5-1) 

Using Eq. ()5-l()|l . it can be rewritten as: 

Wn{x,y) = det ( Id„ + J„(x, ?/)7r„_i^oo0oo*vrJ,-i + vr„,-i^— n„_i— ^— 7r^_i J 

" \ X y r J 



(5-2) 

Wn+i{x, y) = det ( Id„+i + Jn{x, ?/)7r„V^oo0oo*7r* + 7r„ T^^n^i T^-^n ] (5-3) 

n+1 \ X — Q y — J 

The difference is thus the same determinant with substracting 1 in the n, n position: 

Wn+i{x,y) - Wn{x,y) 

det n„_i -|- J ri^ oo^f^ oo 

n+i \ X — Q y — 

(5-4) 



For any arbitrary non- vanishing a and a, multiply the matrix inside the determinant 
Eq. ()5-4p on the left by Id„_|_i — tt^ ^;^^" + (whose determinant is -^)^ and on the 

right by Id„+i — ^"f °° vr* + a^r^ (whose determinant is 

(Id„+i -7r„— ha^— )n„_i(Id„+i — 7r„ + a— — ) = n„_i = Id„+i - e„e„ 

(5-5) 

(Id„+i-7r„^ + a^)J„^o.0oo*(Idn+i-^^7r* +a^) = J„«a^^^^ (5-6) 
(Id„+i - Tin— \- a— — )7r„— — — n„_i 

Yn Yn X 

= (n„_i(l - ^) + a^){R + V^oo0 ooVn-l 

= n„_i(i - ^^)n„i?n„_i + + ae^^^o n„_i 
= n„_i(i - ^:^)n„i?n,_i + ae„T^ 



(5-7) 



where 



Ti = ^ + (fJUn-i , T„ = ^ + n„_i^;o (5-8) 

Vn (Pn 
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Finaly, using Eq. ()5-6|) we have: 



aa 



= det(Id„+i + (J„aa - l)e„e^+ 

n+l 

+ ((n„_i(x - Q)n„_i)"i + ae„T^)((n„_i(|/ - P*)nn-i)"' + «T„e^)) 

(5-9) 



If we choose aa = we have: 

— ^— - - Wn{x, y)) 

det fld„+i + ((n„_i(x - Q)n„„i)-i + ae„T^)((n„_i(y - p*)n„_i)-i 

det fid„ + ((n„_i(y - p*)n„_i)-i + «fX)((n„_i(x - Q)n„_i)-i 
det f id„ + (n„„i(y - p*)n„_i)-i(n„_i(x - Q)Yin-i)'^ + ^t„t^) 

(5-10) 

Then , we multiply both sides by Eq. ()5-l|) . i.e. by det(n„_i(a; — Q)n„_i) and 
det(n„_i(2/-P*)n„_i): 



det (id^ + n„_i(i/ - p*)n„_i(x - g)n„_i 

(5-11) 



+^n„_i(y - p')fnTl,{x - Q)n„_: 



We have: 



T*(x-g)n„_i = ^e*,i?(a;-Q)n„„i + 0jn„_i(x-g)n„_i 

— ej,(i - ^)n„_i + J(x - g) - Jn"(x - g)n„_i 

_4 + 4_0jn"(a;-g)n„_i 
n"(x - g)n„_i 

(5 - 12) 

Notice also that Heisenberg's relation Eq. ()2-45p implies: 

Id„ + n,_i(y - P*)(x - Q)Tln-i = n„_i(x - Q)iy - P*)n„_i (5-13) 
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Thus: 



J n. 



{Wn+i{x,y) - Wn{x,y)) 



det (id, + n„_i(i/ - p*)n„„i(x - Q)n„_i 

n V 

+^n„_i(y - p*)n>;o0oo n"(x - Q)n„_i 
det fid„ + n„_i(i/ - p*)(x - g)n„_i - n„_i(i/ - p*)n"(x - g)n„_i 

n V 



- t 



= det (n„_i(a; - Q)(y - p*)n„_i - n„_i(y - p*)n"(i - ^^^2£f2^)n"(x - Q)n„_i) 

(5-14) 

Now, use Eq. (|5-4p and Eq. (j5-8p . that gives: 

^ {Wn+i{x,y) -Wn{x,y)) 



det (id„ - n„_i(y - p*)n"(i - ^:22^)n"(x - Q)n„_iP*x 



t 



det (id„+, - n„(i - %^)n.pn„_i(y - p*)n"(i - 1^^^ 



(5 - 15) 

Now, because of Eq. ()2-47p . we have: 

n„pn„_i(j/-p*)n" = iinRVl{Q)ir 

" x-Q 
-W„(x) 

(5 - 16) 

Therefore: 



t 



Wn+i{x,y) - Wn{x,y) = det fld„+i - n„(l - ^^^^)W„(1 - ^^^^)W* 

(5-17) 

The second part of the formula is obtained by applying Lemma lL2l 



47 



Appendix F Proof of theorem 18.1 



We start from theorem 17.11 

Wn{x,y) = -fl det {Id n+i- Unix, y)Unix,yy) 

"J ~ ^ ^ 'U 

Let us compute the various terms: 

1 1^ 



i-UnUi 



1 



1 - f/nf/r " ^ 



0{y 

where Z*^'''' is the vector (1, 0, ... , 0)* of dimension ^2 + 1. 
At large y, this expression behaves hke e~^^*^^-* 0{y"'^'^'^^^). 



1 - U^^ri 

" "(t/ri-f/„)A„vi/„ " 

y n n 

y n n 

y -^J{x)^'^{x) 
It behaves as e^^^^^ 0{y^'"'~'^'^~^) at large y. 
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/■(o)*$* I ^ m 

f(o)t \ ^ f(0) 

y n n 

(00) 



Similarly, we have: 



y - Hn{x, x) 



I / 1 \ (00) 



(6-4) 



1-UnUi \x-Hn{y,y) 



--UniJoo = ; ' (6-5) 



but unfortunately, that formula is not sufficient to find the large y behaviour. 
Instead, we prove the following lemma: 



Lemma F.l 



1 - T7 

{x + V^{y))cl>J ^ _ ^ ^ U^^n = -1 - + 0{y-'^-') (6-6) 



and symetricaly: 



{y + Vi{x))cfJiJi ] = -1 - + 0{x-'^-') (6-7) 



proof: 

Compute the large y behaviour of the following matrix: 

W-\iji+{x + Vi)iJn4>J)^n 
= W-\R\X-Q) + {X + Vi)i,n<t>jnn 

= n"-l((^* + ^f^<i)J) {X-Q)+ V'n0oo*(V^2 (y) - Vi{P')) + iJn4>J{Vi{P') + Q))Iin 

= n"-i(^-^(x-Q)+v^„0oo*(^2(^*) + Q))nn (6-8) 

It is easy to see that this expression behaves like 0(|/^^) at large y. The leading term 
in y is: 

-n"-^(x - g)n„ + 4_ie-^^(^)|-(e^^(^)0„(y))e„_ie^ + 0{y-^) 

y oy 
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y Vhn 



y \ In 

(6-9) 



and thus, using C/„ = — ^e„e^_j^ + 0(1), we have: 



(6 - 10) 

and the following determinant is: 



[/„([/* + (x + y2')40oo*)n„ 

-— eX-i ( (x - Q)n„ + — e„_ien + 0{y'^) 

In \ In J 

In \ In / 



It follows: 



det(l - Un (Ul + ix + V^)i^r^<^>oo')Tln) 

^ det(l + -e, (ei_,{x - Q)n^ + - e^) + 0{y-')) 

7n V 7n / 

= l + -fe*_i(x-g)e„ + -)+0(y-^)) 

. i + (i-i)^0(.-)) 

= ^ + o(y-')) (6-11) 



4 + 0(1/-^)) 

det(l - (t/* + {X + ^')^n0oo*)n 



det(l -[/„[/* + (X + ^')^n ^„0oo*)nn) 

det(l - [/„ L^*) det(l + (x + V^)- l-^Un V'n0oo*)n„) 

det(l-t/„i/*) fl + (a; + V,O0oo \ -. ^nV^n) 



(6 - 12) 
that implies: 

1 ^ 77/ 
(x + y^O'/'oo*— — — -f/nV'n = -1 + — -— -^+0(y-''^-') 



1 7^det(l-t/„f/*) 

n 

^n(a^,?/) 
n 
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-l-7r7% + 0(t/-'^^-^) 



-1 - + ©(y-'^-^) 



(6 - 13) 



similarly we have: 

(y + V;{x))<fjui I ^n = -l- + 0{x-'^-') (6-14) 

□ 



F.5 U-conjecture. Lemma 18.11 
Using lemma IFTTI we have: 



(x + Vi{y)) det(l - [/„ ill) {cfJiJi _ ^ ij^) (0OO* _ I ~ Un ^n) 

1 U Yi U ^ 1 Ufi u ^ 

- t ~ 1 72 - t ~ 1 

- det(l - f/„ Ui) Ui _ ^„) - - Ui _ i^n) 

~, ( 1 y°°) n ( 1 y°°^ 

det(l f/nf/„) (^^ _ ^^(^3,^3,) J ^2 (^y _ if„(a;,x)y 

+o(r') 

i / 1 n ( 1 



— det(y-if„(x,a;)) . „ ^ ^ 1 2 \ u t ) 

+0{y-') 



t / 1 



\ (00) 

-det(y-i/„(x,x)) ( -] +0{y-') 

In \y - Hn{x,x)J 

(6 - 15) 

The first term is a polynomial in y, and this proves theorem 18.11 



Appendix G Proof of theorem 18.3 

Starting from lemma lF\T| we have: 



det(l - f/„f/*) {y + Vlix)) (<fjui I 

Tl 

= ^n{x,y){l + ^^^) + 0{x-') 
Sn{x,y) + n + 0{x-') 

(7-1) 



from which theorem 18.31 follows. 
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Appendix H Proof of theorem 18.2 



Let r„(x) be the matrix: 

1 

7n+l 

By definition, it is such that: 



r„(x) = K-i-+'F^{x) = U:-t-+' + e„+ie^(x - Q) (8-1) 



Thus: 



=r„(x)^„(x) (8-2) 



I?„+i(x) = K+iKli = (r„^„)'(r„^„)-^ 
(8-3) 



Therefore: 



£rr+i = idetiy-V^-r-X) (8-4) 
Now, notice that r~^r' = —^r-^ — e„_d„e* , thus 

Sn+i -£n = {-If-^^ det(C„) (8-5) 

'^n,n—d2 

where we have defined 

C^{x,y) =Ul-''-+\y + V^{x))Kzf- (8-6) 

We have: 

Cn = K''''^\y - P')K~-i' + nr'^^+i^^M^^^(i - n„_i)(x - Q)Kzt' 

(8-7) 

and: 

CnK-t'R'K'^'^^ = nr^2+i(i - w„w*)nr'''+' (8-8) 

That imphes: 

det(C„) UlztR'K'''^' = (-l)''7n . . . 7n-d,+i det(l - W„W* ) 



(8-9) 

which proves the formula. 



^y,+iQn^ det(l-W„>V*; 
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Appendix I A usefull formula for determinants with 
rank 2 matrices 

Lemma 1.2 if M is an invertible matrix, and a, b, c, d are arbitrary vectors, we have: 



det(M + a6* + cd*) = detM { (1 + b'—a)(l + d'—c) - b'—c d'—a ] (9-1) 
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